Uncovering the mesendoderm gene regulatory network through multi-omic data integration.

Mesendodermal specification is one of the earliest events in embryogenesis, where cells first acquire distinct identities. Cell differentiation is a highly regulated process that involves the function of numerous transcription factors (TFs) and signaling molecules, which can be described with gene regulatory networks (GRNs). Cell differentiation GRNs are difficult to build because existing mechanistic methods are low throughput, and high-throughput methods tend to be non-mechanistic. Additionally, integrating highly dimensional data composed of more than two data types is challenging. Here, we use linked self-organizing maps to combine chromatin immunoprecipitation sequencing (ChIP-seq)/ATAC-seq with temporal, spatial, and perturbation RNA sequencing (RNA-seq) data from Xenopus tropicalis mesendoderm development to build a high-resolution genome scale mechanistic GRN. We recover both known and previously unsuspected TF-DNA/TF-TF interactions validated through reporter assays. Our analysis provides insights into transcriptional regulation of early cell fate decisions and provides a general approach to building GRNs using highly dimensional multi-omic datasets.

Tbx5 drives Aldh1a2 expression to regulate a RA-Hedgehog-Wnt gene regulatory network coordinating cardiopulmonary development.

The gene regulatory networks that coordinate the development of the cardiac and pulmonary systems are essential for terrestrial life but poorly understood. The T-box transcription factor Tbx5 is critical for both pulmonary specification and heart development, but how these activities are mechanistically integrated remains unclear. Here using Xenopus and mouse embryos, we establish molecular links between Tbx5 and retinoic acid (RA) signaling in the mesoderm and between RA signaling and sonic hedgehog expression in the endoderm to unveil a conserved RA-Hedgehog-Wnt signaling cascade coordinating cardiopulmonary (CP) development. We demonstrate that Tbx5 directly maintains expression of aldh1a2, the RA-synthesizing enzyme, in the foregut lateral plate mesoderm via an evolutionarily conserved intronic enhancer. Tbx5 promotes posterior second heart field identity in a positive feedback loop with RA, antagonizing a Fgf8-Cyp regulatory module to restrict FGF activity to the anterior. We find that Tbx5/Aldh1a2-dependent RA signaling directly activates shh transcription in the adjacent foregut endoderm through a conserved MACS1 enhancer. Hedgehog signaling coordinates with Tbx5 in the mesoderm to activate expression of wnt2/2b, which induces pulmonary fate in the foregut endoderm. These results provide mechanistic insight into the interrelationship between heart and lung development informing CP evolution and birth defects.

The dorsal blastopore lip is a source of signals inducing planar cell polarity in the Xenopus neural plate.

Coordinated polarization of cells in the tissue plane, known as planar cell polarity (PCP), is associated with a signaling pathway critical for the control of morphogenetic processes. Although the segregation of PCP components to opposite cell borders is believed to play a critical role in this pathway, whether PCP derives from egg polarity or preexistent long-range gradient, or forms in response to a localized cue, remains a challenging question. Here we investigate the Xenopus neural plate, a tissue that has been previously shown to exhibit PCP. By imaging Vangl2 and Prickle3, we show that PCP is progressively acquired in the neural plate and requires a signal from the posterior region of the embryo. Tissue transplantations indicated that PCP is triggered in the neural plate by a planar cue from the dorsal blastopore lip. The PCP cue did not depend on the orientation of the graft and was distinct from neural inducers. These observations suggest that neuroectodermal PCP is not instructed by a preexisting molecular gradient but induced by a signal from the dorsal blastopore lip.

Spatio-temporal mRNA tracking in the early zebrafish embryo.

Early stages of embryogenesis depend on subcellular localization and transport of maternal mRNA. However, systematic analysis of these processes is hindered by a lack of spatio-temporal information in single-cell RNA sequencing. Here, we combine spatially-resolved transcriptomics and single-cell RNA labeling to perform a spatio-temporal analysis of the transcriptome during early zebrafish development. We measure spatial localization of mRNA molecules within the one-cell stage embryo, which allows us to identify a class of mRNAs that are specifically localized at an extraembryonic position, the vegetal pole. Furthermore, we establish a method for high-throughput single-cell RNA labeling in early zebrafish embryos, which enables us to follow the fate of individual maternal transcripts until gastrulation. This approach reveals that many localized transcripts are specifically transported to the primordial germ cells. Finally, we acquire spatial transcriptomes of two xenopus species and compare evolutionary conservation of localized genes as well as enriched sequence motifs.

The dual-specificity protein kinase Clk3 is essential for Xenopus neural development.

During vertebrate development, the formation of the central nervous system (CNS) is initiated by neural induction and patterning of the embryonic ectoderm. We previously reported that Cdc2-like kinase 2 (Clk2) promotes neural development in Xenopus embryos by regulating morphogen signaling. However, the functions of other Clk family members and their roles in early embryonic development remain unknown. Here, we show that in addition to Clk2, Clk1 and Clk3 play a role in the formation of neural tissue in Xenopus. clk1 and clk3 are co-expressed in the developing neural tissue during early Xenopus embryogenesis. We found that overexpression of clk1 and clk3 increases the expression of neural marker genes in ectodermal explants. Furthermore, knockdown experiments showed that clk3 is required for the formation of neural tissues. These results suggest that Xenopus Clk3 plays an essential role in promoting neural development during early embryogenesis.

Mechanical stretch scales centriole number to apical area via Piezo1 in multiciliated cells.

How cells count and regulate organelle number is a fundamental question in cell biology. For example, most cells restrict centrioles to two in number and assemble one cilium; however, multiciliated cells (MCCs) synthesize hundreds of centrioles to assemble multiple cilia. Aberration in centriole/cilia number impairs MCC function and can lead to pathological outcomes. Yet how MCCs control centriole number remains unknown. Using Xenopus, we demonstrate that centriole number scales with apical area over a remarkable 40-fold change in size. We find that tensile forces that shape the apical area also trigger centriole amplification based on both cell stretching experiments and disruption of embryonic elongation. Unexpectedly, Piezo1, a mechanosensitive ion channel, localizes near each centriole suggesting a potential role in centriole amplification. Indeed, depletion of Piezo1 affects centriole amplification and disrupts its correlation with the apical area in a tension-dependent manner. Thus, mechanical forces calibrate cilia/centriole number to the MCC apical area via Piezo1. Our results provide new perspectives to study organelle number control essential for optimal cell function.

Cell cycle control during early embryogenesis.

Understanding the mechanisms of embryonic cell cycles is a central goal of developmental biology, as the regulation of the cell cycle must be closely coordinated with other events during early embryogenesis. Quantitative imaging approaches have recently begun to reveal how the cell cycle oscillator is controlled in space and time, and how it is integrated with mechanical signals to drive morphogenesis. Here, we discuss how the Drosophila embryo has served as an excellent model for addressing the molecular and physical mechanisms of embryonic cell cycles, with comparisons to other model systems to highlight conserved and species-specific mechanisms. We describe how the rapid cleavage divisions characteristic of most metazoan embryos require chemical waves and cytoplasmic flows to coordinate morphogenesis across the large expanse of the embryo. We also outline how, in the late cleavage divisions, the cell cycle is inter-regulated with the activation of gene expression to ensure a reliable maternal-to-zygotic transition. Finally, we discuss how precise transcriptional regulation of the timing of mitosis ensures that tissue morphogenesis and cell proliferation are tightly controlled during gastrulation.

Three-dimensional folding dynamics of the Xenopus tropicalis genome.

Animal interphase chromosomes are organized into topologically associating domains (TADs). How TADs are formed is not fully understood. Here, we combined high-throughput chromosome conformation capture and gene silencing to obtain insights into TAD dynamics in Xenopus tropicalis embryos. First, TAD establishment in X. tropicalis is similar to that in mice and flies and does not depend on zygotic genome transcriptional activation. This process is followed by further refinements in active and repressive chromatin compartments and the appearance of loops and stripes. Second, within TADs, higher self-interaction frequencies at one end of the boundary are associated with higher DNA occupancy of the architectural proteins CTCF and Rad21. Third, the chromatin remodeling factor ISWI is required for de novo TAD formation. Finally, TAD structures are variable in different tissues. Our work shows that X. tropicalis is a powerful model for chromosome architecture analysis and suggests that chromatin remodeling plays an essential role in de novo TAD establishment.

Mechanical heterogeneity along single cell-cell junctions is driven by lateral clustering of cadherins during vertebrate axis elongation.

Morphogenesis is governed by the interplay of molecular signals and mechanical forces across multiple length scales. The last decade has seen tremendous advances in our understanding of the dynamics of protein localization and turnover at subcellular length scales, and at the other end of the spectrum, of mechanics at tissue-level length scales. Integrating the two remains a challenge, however, because we lack a detailed understanding of the subcellular patterns of mechanical properties of cells within tissues. Here, in the context of the elongating body axis of Xenopus embryos, we combine tools from cell biology and physics to demonstrate that individual cell-cell junctions display finely-patterned local mechanical heterogeneity along their length. We show that such local mechanical patterning is essential for the cell movements of convergent extension and is imparted by locally patterned clustering of a classical cadherin. Finally, the patterning of cadherins and thus local mechanics along cell-cell junctions are controlled by Planar Cell Polarity signaling, a key genetic module for CE that is mutated in diverse human birth defects.

Combinatorial transcription factor activities on open chromatin induce embryonic heterogeneity in vertebrates.

During vertebrate gastrulation, mesoderm is induced in pluripotent cells, concomitant with dorsal-ventral patterning and establishing of the dorsal axis. We applied single-cell chromatin accessibility and transcriptome analyses to explore the emergence of cellular heterogeneity during gastrulation in Xenopus tropicalis. Transcriptionally inactive lineage-restricted genes exhibit relatively open chromatin in animal caps, whereas chromatin accessibility in dorsal marginal zone cells more closely reflects transcriptional activity. We characterized single-cell trajectories and identified head and trunk organizer cell clusters in early gastrulae. By integrating chromatin accessibility and transcriptome data, we inferred the activity of transcription factors in single-cell clusters and tested the activity of organizer-expressed transcription factors in animal caps, alone or in combination. The expression profile induced by a combination of Foxb1 and Eomes most closely resembles that observed in the head organizer. Genes induced by Eomes, Otx2, or the Irx3-Otx2 combination are enriched for maternally regulated H3K4me3 modifications, whereas Lhx8-induced genes are marked more frequently by zygotically controlled H3K4me3. Taken together, our results show that transcription factors cooperate in a combinatorial fashion in generally open chromatin to orchestrate zygotic gene expression.

The molecular dynamics of subdistal appendages in multi-ciliated cells.

The motile cilia of ependymal cells coordinate their beats to facilitate a forceful and directed flow of cerebrospinal fluid (CSF). Each cilium originates from a basal body with a basal foot protruding from one side. A uniform alignment of these basal feet is crucial for the coordination of ciliary beating. The process by which the basal foot originates from subdistal appendages of the basal body, however, is unresolved. Here, we show FGFR1 Oncogene Partner (FOP) is a useful marker for delineating the transformation of a circular, unpolarized subdistal appendage into a polarized structure with a basal foot. Ankyrin repeat and SAM domain-containing protein 1A (ANKS1A) interacts with FOP to assemble region I of the basal foot. Importantly, disruption of ANKS1A reduces the size of region I. This produces an unstable basal foot, which disrupts rotational polarity and the coordinated beating of cilia in young adult mice. ANKS1A deficiency also leads to severe degeneration of the basal foot in aged mice and the detachment of cilia from their basal bodies. This role of ANKS1A in the polarization of the basal foot is evolutionarily conserved in vertebrates. Thus, ANKS1A regulates FOP to build and maintain the polarity of subdistal appendages.

Nutrient availability contributes to a graded refractory period for regeneration in Xenopus tropicalis.

Xenopus tadpoles are a unique model for regeneration in that they exhibit two distinct phases of age-specific regenerative competence. In Xenopus laevis, young tadpoles fully regenerate following major injuries such as tail transection, then transiently lose regenerative competence during the "refractory period" from stages 45-47. Regenerative competence is then regained in older tadpoles before being permanently lost during metamorphosis. Here we show that a similar refractory period exists in X. tropicalis. Notably, tadpoles lose regenerative competence gradually in X. tropicalis, with full regenerative competence lost at stage 47. We find that the refractory period coincides closely with depletion of maternal yolk stores and the onset of independent feeding, and so we hypothesized that it might be caused in part by nutrient stress. In support of this hypothesis, we find that cell proliferation declines throughout the tail as the refractory period approaches. When we block nutrient mobilization by inhibiting mTOR signaling, we find that tadpole growth and regeneration are reduced, while yolk stores persist. Finally, we are able to restore regenerative competence and cell proliferation during the refractory period by abundantly feeding tadpoles. Our study argues that nutrient stress contributes to lack of regenerative competence and introduces the X. tropicalis refractory period as a valuable new model for interrogating how metabolic constraints inform regeneration.

Nucleoporin NUP205 plays a critical role in cilia and congenital disease.

Ciliopathies affect a variety of tissues during development including the heart, kidneys, respiratory tract, and retina. Though an increasing number of monogenic causes of ciliopathies have been described, many remain unexplained. Recently, recessive variants in NUP93 and NUP205 encoding two proteins of the inner ring of the nuclear pore complex were implicated as causes of steroid resistant nephrotic syndrome. In addition, we previously found that the inner ring nucleoporins NUP93 and NUP188 function in proper left-right patterning in developing embryos via a role at the cilium. Here, we describe the role of an additional inner ring nucleoporin NUP205 in cilia biology and establishment of normal organ situs. Using knockdown in Xenopus, we show that Nup205 depletion results in loss of cilia and abnormal cardiac morphology. Furthermore, by transmission electron microscopy, we observe a loss of cilia and mispositioning of intracellular ciliary structures such as basal bodies and rootlets upon depleting inner ring nucleoporins. We describe a model wherein NUP93 interacting with either NUP188 or NUP205 is necessary for cilia. We thus provide evidence that dysregulation of inner ring nucleoporin genes that have been identified in patients may contribute to pathogenesis through cilia dysfunction.

Analysis of Thyroid Hormone Receptor α-Knockout Tadpoles Reveals That the Activation of Cell Cycle Program Is Involved in Thyroid Hormone-Induced Larval Epithelial Cell Death and Adult Intestinal Stem Cell Development During Xenopus tropicalis Metamorphosis.

Background: There are two highly conserved thyroid hormone (triiodothyronine [T3]) receptor (TR) genes, TRα and TRß, in all vertebrates, and the expression of TRα but not TRß is activated earlier than T3 synthesis during development. In human, high levels of T3 are present during the several months around birth, and T3 deficiency during this period causes severe developmental abnormalities including skeletal and intestinal defects. It is, however, difficult to study this period in mammals as the embryos and neonates depend on maternal supply of nutrients for survival. However, Xenopus tropicalis undergoes a T3-dependent metamorphosis, which drastically changes essentially every organ in a tadpole. Of interest is intestinal remodeling, which involves near complete degeneration of the larval epithelium through apoptosis. Concurrently, adult intestinal stem cells are formed de novo and subsequently give rise to the self-renewing adult epithelial system, resembling intestinal maturation around birth in mammals. We have previously demonstrated that T3 signaling is essential for the formation of adult intestinal stem cells during metamorphosis. Methods: We studied the function of endogenous TRα in the tadpole intestine by using knockout animals and RNA-seq analysis. Results: We observed that removing endogenous TRα caused defects in intestinal remodeling, including drastically reduced larval epithelial cell death and adult intestinal stem cell proliferation. Using RNA-seq on intestinal RNA from premetamorphic wild-type and TRα-knockout tadpoles treated with or without T3 for one day, before any detectable T3-induced cell death and stem cell formation in the tadpole intestine, we identified more than 1500 genes, which were regulated by T3 treatment of the wild-type but not TRα-knockout tadpoles. Gene Ontology and biological pathway analyses revealed that surprisingly, these TRα-regulated genes were highly enriched with cell cycle-related genes, in addition to genes related to stem cells and apoptosis. Conclusions: Our findings suggest that TRα-mediated T3 activation of the cell cycle program is involved in larval epithelial cell death and adult epithelial stem cell development during intestinal remodeling.

The RNA helicase DDX3 induces neural crest by promoting AKT activity.

Mutations in the RNA helicase DDX3 have emerged as a frequent cause of intellectual disability in humans. Because many individuals carrying DDX3 mutations have additional defects in craniofacial structures and other tissues containing neural crest (NC)-derived cells, we hypothesized that DDX3 is also important for NC development. Using Xenopus tropicalis as a model, we show that DDX3 is required for normal NC induction and craniofacial morphogenesis by regulating AKT kinase activity. Depletion of DDX3 decreases AKT activity and AKT-dependent inhibitory phosphorylation of GSK3ß, leading to reduced levels of ß-catenin and Snai1: two GSK3ß substrates that are crucial for NC induction. DDX3 function in regulating these downstream signaling events during NC induction is likely mediated by RAC1, a small GTPase whose translation depends on the RNA helicase activity of DDX3. These results suggest an evolutionarily conserved role of DDX3 in NC development by promoting AKT activity, and provide a potential mechanism for the NC-related birth defects displayed by individuals harboring mutations in DDX3 and its downstream effectors in this signaling cascade.

Hes5.9 Coordinate FGF and Notch Signaling to Modulate Gastrulation via Regulating Cell Fate Specification and Cell Migration in Xenopus tropicalis.

Gastrulation drives the establishment of three germ layers and embryonic axes during frog embryonic development. Mesodermal cell fate specification and morphogenetic movements are vital factors coordinating gastrulation, which are regulated by numerous signaling pathways, such as the Wnt (Wingless/Integrated), Notch, and FGF (Fibroblast growth factor) pathways. However, the coordination of the Notch and FGF signaling pathways during gastrulation remains unclear. We identified a novel helix-loop-helix DNA binding domain gene (Hes5.9), which was regulated by the FGF and Notch signaling pathways during gastrulation. Furthermore, gain- and loss-of-function of Hes5.9 led to defective cell migration and disturbed the expression patterns of mesodermal and endodermal marker genes, thus interfering with gastrulation. Collectively, these results suggest that Hes5.9 plays a crucial role in cell fate decisions and cell migration during gastrulation, which is modulated by the FGF and Notch signaling pathways.

A nonsense variant in Rap Guanine Nucleotide Exchange Factor 5 (RAPGEF5) is associated with equine familial isolated hypoparathyroidism in Thoroughbred foals.

Idiopathic hypocalcemia in Thoroughbred (TB) foals causes tetany and seizures and is invariably fatal. Based upon the similarity of this disease with human familial hypoparathyroidism and occurrence only in the TB breed, we conducted a genetic investigation on two affected TB foals. Familial hypoparathyroidism was identified, and pedigree analysis suggested an autosomal recessive (AR) mode of inheritance. We performed whole-genome sequencing of the two foals, their unaffected dams and four unaffected, unrelated TB horses. Both homozygosity mapping and an association analysis were used to prioritize potential genetic variants. Of the 2,808 variants that significantly associated with the phenotype using an AR mode of inheritance (P<0.02) and located within a region of homozygosity, 1,507 (54%) were located in a 9.7 Mb region on chr4 (44.9-54.6 Mb). Within this region, a nonsense variant (RAPGEF5 c.2624C>A,p.Ser875*) was significantly associated with the hypoparathyroid phenotype (Pallelic = 0.008). Affected foals were homozygous for the variant, with two additional affected foals subsequently confirmed in 2019. Necropsies of all affected foals failed to identify any histologically normal parathyroid glands. Because the nonsense mutation in RAPGEF5 was near the C-terminal end of the protein, the impact on protein function was unclear. Therefore, we tested the variant in our Xenopus overexpression model and demonstrated RAPGEF5 loss-of-function. This RAPGEF5 variant represents the first genetic variant for hypoparathyroidism identified in any domestic animal species.

Thyroid Hormone Induces DNA Demethylation in Xenopus Tadpole Brain.

Thyroid hormone (T3) plays pivotal roles in vertebrate development, acting via nuclear T3 receptors (TRs) that regulate gene transcription by promoting post-translational modifications to histones. Methylation of cytosine residues in deoxyribonucleic acid (DNA) also modulates gene transcription, and our recent finding of predominant DNA demethylation in the brain of Xenopus tadpoles at metamorphosis, a T3-dependent developmental process, caused us to hypothesize that T3 induces these changes in vivo. Treatment of premetamorphic tadpoles with T3 for 24 or 48 hours increased immunoreactivity in several brain regions for the DNA demethylation intermediates 5-hydroxymethylcytosine (5-hmC) and 5-carboxylcytosine, and the methylcytosine dioxygenase ten-eleven translocation 3 (TET3). Thyroid hormone treatment induced locus-specific DNA demethylation in proximity to known T3 response elements within the DNA methyltransferase 3a and Krüppel-like factor 9 genes, analyzed by 5-hmC immunoprecipitation and methylation sensitive restriction enzyme digest. Chromatin-immunoprecipitation (ChIP) assay showed that T3 induced TET3 recruitment to these loci. Furthermore, the messenger ribonucleic acid for several genes encoding DNA demethylation enzymes were induced by T3 in a time-dependent manner in tadpole brain. A TR ChIP-sequencing experiment identified putative TR binding sites at several of these genes, and we provide multiple lines of evidence to support that tet2 contains a bona fide T3 response element. Our findings show that T3 can promote DNA demethylation in developing tadpole brain, in part by promoting TET3 recruitment to discrete genomic regions, and by inducing genes that encode DNA demethylation enzymes.